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ABSTRACT 

Precision cosmology studies based on wide-field surveys of galaxy clusters 
benefit from constraints on intrinsic scatter in mass-observable relationships. In 
principle, two-parameter models combining direct measurements of galaxy cluster 
structural variation with mass proxies such as X-ray luminosity and temperature 
can be used to constrain scatter in the relationship between the mass proxy and 
the cluster's halo mass and to measure the redshift evolution of that scatter. One 
candidate for quantifying cluster substructure is the ICM temperature inhomo- 
geneity inferred from X-ray spectral properties, an example of which is Thbr, 
the ratio of hardband to broadband spectral-fit temperatures. In this paper we 
test the effectiveness of Thbr as an indicator of scatter in the mass-temperature 
relation using 118 galaxy clusters simulated with radiative cooling and feedback. 
We find that, while T H br is correlated with clusters' departures 5\nT x from the 
mean M-Tx relation, the effect is modest. 

Subject headings: galaxies: clusters: general, X-rays: galaxies: clusters 



1. Introduction 



Galaxy clusters play an important role in precision cosmology that com plements other 
techniques like Type la supernovae luminosity-distance relation measurements (IPerlmutter et al 
19991 : iRiess et al.lll998l). baryonic acoustic oscillations (BAOs) angular-distance relation mea- 



surements ( 
diation (see 



Frieman et al. 



isenstein et al. 120051) . and observations of the cos mic microwave backgro und ra 



( 120081 ) for a review). For example, IVikhlinin et al.l (120091 ) exploit 
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dark energy's influence on the growth of structure by using 37 moderate-redshift and 49 low- 
redshift clusters to measure the shape of the galaxy-cluster mass function and its redshift 
evolution, which constrain the dark-energy density parameter Q\ to 0.83 ±0.15 in a non-flat 
ACDM cosmol ogy and the dark en ergy equation of state parameter wq to —1.14 ± 0.21 in a 
flat cosmology. iMantz et al.l (120101 ) have obtained similar results from measurements of the 
evolving number density of the largest clusters in order to constrain wq to —1.01 ± 0.20. 

Strategies like these that compare model predictions to galaxy-cluster sample statistics 
inevitably confront sample error and sample bias. Future surveys expected to gather samples 
of 10-40 thousand galaxy cluster^ will maximize survey reach while maintaining sufficient 
observation quality in order to minimize sample error, but they still must grapple with a 
major source of sample bias, which is scatter in the relationship used to infer cluster mass 
from an observable mass proxy. An important mass-observable relation for galaxy cluster 
studies connects dark matter halo mass to the temperature of the intracluster medium (ICM) 
inferred from its X-ray spectrum (its "X-ray temperature" Tx)- In this paper we investigate 
the possibility of correcting for scatter in this relation using temperature-inhomogeneity in 
the ICM, and discuss challenges that may exist in such a program. 

A significant amount of uncertainty in the dark-energy constraints obtainable from large 
cluster surveys derives from uncertai nty in scatter about the mean scaling relat ions obeyed 
by galaxy clusters' bulk properties ( ILima &: Hul l2005t ICunha &: Evrardl I201CH ) . The key 
galaxy cluster property to measure when trying to constrain dark energy with clusters is 
the cluster's mass, which cannot be directly observed. Theoretical considerations predict 
correlations among halo mass and more readily observed cluster properties, like its galaxy 
richness, the velocity dispersion of its galaxies, Tx, the Sunyaev-Zel'dovich decrement, the 
gas mass, and Yx p arameter, which is th e product of Tx and the gas mass inferred from 
X-ray observations ( jKravtsov et al.ll2006l ). Theory also predicts intrin sic scatter in these 



Fab j an et al. 


2011; 


Rasia et al. 


2011) 



One way to deal with intrinsic scatter is to join a cluster model to a cosmological model 
and simultaneously fit for the parameters of both, leveraging the statistical power of large 
surveys and "self-calibrating" the mass-observable relations. Another approach is to combine 
observables that tend to depart from the expected scaling relations in opposite ways, yielding 
a new, low-scatter composite observable. An example low-scatter composite observable is 
lx, since at a given halo mass, offsets in the measured g as mass at fixed tota l mass tend to 
anti-correlate with offsets in the measured temperature (jKravtsov et al.l 120061 ) . 



1 http://www. mpe.mpg.de/heg/www/Projects/EROSITA/main. html 
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Another family of low-scatter composite observables attempt to measure structural vari- 
ation directly. Mergers, relaxation, and non-adiabatic processes like radiative cooling, star 
formation, and feedback ought to leave a visible imprint that may allow us to measure and 
correct for sc atter. For example, o ne might use imaging to quantify resolved morphological 
substructure. IJeltema et al.l (120081) ap ply two observationally - motivated struct ure measures, 



the power ratios (IBuote fc Tsailll995l ) an d the centroid shift (IMohr et al.lll993l), to a sample 



of galaxy clusters simulated with Enzo (INorman fc Bryan! 11999 



O'Shea et al. 



200J). They 



find that cluster structure correlates strongly with bias in mass estimates derived from Tx 
under the assumption of hydrostatic equili brium and account i ng for cluster structure can be 
used to correct some of the bi as. Similarly, Ventimiglia et al.l ( 120081 ) apply the power ratios, 



centroid shift, and axial ratio (jO'Hara et al 



20061 ) substructure measures to the same sample 



of simulated clusters used in this paper, and find that cluster substructure correlates with de- 
partures from the mean M-Tx relationship in the sense that clusters with more substructure 
tend to have a lower temperature at a give n mass, and can be used to refi ne mass estimates 
derived from the ICM X-ray temperature. iPiffaretti fc Valdarninil (120081 ) likewise find that 
greater subs tructure, as quanti fied with power ratios, correlates with lower temperature at a 



given mass. lYang et al.l (120091 ) find a strong correlation between mass-temperature scatter 
and halo concentration in their sample of simulated clusters, with cooler clusters appearing 
more concentrated than warmer clusters at similar mass. 



Jeltema et al.l ( 120081 ) observe, however, that line-of-sight projection effects lead to sig- 
nificant uncertainties in morphologically-derived substructure measures. Substructure also 
becomes more difficult to resolve at high redshift. Spectral signatures of dynamical state are 
therefore attractive because they are aspect-independent and redshift-in dependent. One such 
spect r al signature of dynami cal state is the "temperature ratio" Thbr ( jMathiesen fc Evrard 
200 ll ; ICavagnolo et al.l [20081) , which divides a "hardband" spectral-fit temperature by a 
"broadband" spectral-fit temperature. An energy cut applied to a broadband spectrum 
produces a hardband spectrum and s erves to filter out cooler line-emitting components of 
the ICM. iMathiesen fc Evrardl (j200lf ) studied the effects of relaxation on the observable 
properties of galaxy clusters, usi ng a sample of numerically-simulated galaxy clusters gen- 
erated by IMohr fc Evrardl ( 119971 ). They found that hardband (2.0-9.0 keV) X-ray spectral 
fit temperatures average ~20% higher than broadband (0.5-9.0 keV) temperatures and sug- 
gested that this effect m ay signal the pres ence of cool, luminous sub-clusters lowering the 
broadband temperature. IValdarninil (120061 ) corroborated these findings in simulations that 
included radiative cooling. 

Figures[TJand|2]illustrate the effect. Both show a counts spectrum and single-temperature 
fit for a typical, unrelaxed cluster in our sample, with an aperture set at i?25oo and a core 
region out to 0.15i?25oo excised. Figured] is for a single-temperature model fit to the broad 
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band, while figure |2] is for the hard band. Note the excess emission relative to the model 
above 4.0 keV and the deficit below 2.5 keV that arises because the model cannot simul- 



taneously fit both a hot component and a cooler, line-emitting component ( IMazzotta et al. 



20041 ) . Figure [2] shows the same spectrum, but with the residuals for a model fit just to 
the hard band. In this case the fit is much better over the range from 2.0 to 7.0 keV but 
under-predicts the emission below 2.0 keV. 



Mathiesen fc Evrardl (120011 ) suggested that the temperature skewing the y observed might 



i ndica te a real temperature skewing detectable in real clusters using Chandra. ICavagnolo et al. 
(120081 ) fit single-temperature emission models to the hard band (2.0-7.0 keV) and broad band 
(0.7-7.0 keV) for a large (N = 192) sample of clusters with observations selected from the 
Chandra Data Archive. These authors show that for a large, heterogeneous sample of clus- 
ters across a broad temperature range, the distribution of Thbr has a mean of 1.16 and 
an rms deviation a = ±0.10, with Thbr tending to be larger in merging systems. They 
also report that while this signal is significant in the aggregate, the errors for any single 
^hbr measurement and the scatter across all of the measurements together pose challenges 
for any effort to use T H br either to select for merging systems or to obtain more accurate 
mass estimates. In order to meet these challenges it is important also to examine Thbr 
and similar spectral signatures of dynam ical state in a simulation context. Like the origi- 
nal study of iMathiesen &: Evrardl (120011 ). this paper examines the temperature ratio for a 
sample of simulated c luster s and simulated Chandra observations, and as in the subsequent 
study by IValdarnini (l2006h . the simulated clusters analyzed here were generated using a 
hydrodynamical code with radiative cooling included. 

This paper is organized as follows. In §2, we describe the simulated clusters in our 
sample along with the X-MAS code for generating their mock X-ray observations. In §3, 
we present our analysis methods and define the spectral signatures of dynamical state we 
examine. In §4, we report and discuss our results, and §5 summarizes our work. 



2. Methods 



2.1. Numerical Simulations 



This study is based on an analysis of 118 clusters simul ated using the cosmological 
hydrodynamics TREE+SPH code GADGET-2 (jSpringell 120051 ) . which were simulated in a 
standard A cold dark matter (ACDM) universe with matter density Qm — 0.3, h = 0.7, fi& 
= 0.04, and cr§ = 0.8. The simulation includes radiative cooling assuming an optically-thin 
gas of primordial composition, with a time- dependent UV background from a population of 
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Fig. 1. — X-MAS simulated counts spectrum for a simulated cluster that appears to have 
significant temperature structure. A single-temperature MEKAL model fit to the [0.7-7.0] 
keV broad band is over-plotted as the solid line, with ksT = 2.42 ± 0.03 keV. Fit residuals 
appear in the bottom panel. This figure and Figure |2] illustrate qualitatively the effect that 
additional cool components have on a single-temperature fit. 
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Fig. 2. — X-MAS simulated counts spectrum for the same simulated cluster as in the previous 
figure. A single-temperature MEKAL model fit to the [2.0-7.0] keV hard band is over-plotted 
as the solid line, with ksT = 3.52 ± 0.2 keV. Fit residuals appear in the bottom panel. 
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quasars, and handles star formation and supernova feedback using a two-phase fluid model 
with cold star-forming clouds em bedded in a hot medium. All but four of the clusters are 
from the simulation described in iBorganil (120041 ). who simulated a box 192 h~ l Mpc on a 
side, with 480 3 dark matter particles and an equal number of gas particles. The present 
analysis considers the 114 most massive clusters within this box at z — 0, which all have 
M200 greater than 5 x 10 13 h~ 1 M & . These are referred to as the B04 Sample in the remainder 
of this paper. By convention, Ma refers to the mass contained in a sphere which has a mean 
density of A times the critical density p c , and whose radius is denoted by 

That cluster set covers the ~1.5-5 keV temperature range, but the 192 h~ x Mpc box is 
too small to contain significantly hotter clusters. We therefore supplemented it with four 
clusters with masses > 10 15 h~ l M & and temp eratures > 5 keV drawn from a dark-matter- 
only simulation in a larger 479 hr 1 Mpc box (IDolag fc Stasyszynl l2009i ). referred to in this 
paper as the D09 sample. The cosmology for this simulation also was ACDM, but with 
erg = 0.9. These were then re-simulated including hydrodynamics, radiative cooling, and 
star formation , again with GADGET-2 and using the zoomed-initial-conditions technique of 



Tormenl (119971 ). with a fourfold increase in resolution. This is comparable to the resolution 



of the clusters in the smaller box. Adding these four massive clusters to our sample gives 
a total of 118 clusters with M 200 in the range 5 x 10 13 h- l M & to 2 x 10 15 h~ l M Q . The 
mean structural properties of massive clusters drawn from a sample with a$ = 0.9 may differ 
somewhat from those of similar- mass clusters in a a 8 = 0.8 universe because they reflect 
a more advanced state of cosmic evolution. However, these four additional clusters carry 
minimal statistical weight in the context of the overall sample. They are included primarily 
to evaluate whether the mean and dispersion of their Thbr values are consistent with those 
of the lower-mass systems. 



2.2. X-Ray Simulations 



The simulated galaxy clusters in our sample are proce ssed with the X-ray Map Simulator 
version 2 (X-MAS) (jGardini et al.ll2004l ; iRasia et al.ll2008l ) to generate X-ray images suitable 
for standard Chandra reduction techniques. In its first step X-MAS uses the outputs of the 
hydrodynamic code to calculate the emissivity of each simulation element within the chosen 
field of view and to project this onto the image plane. In its second step it convolves the 
resulting flux with the appropriate response of a given detector. In the case of our simulated 
Chandra observations, the second step applies the response matrix file and ancillary response 
file for the ACIS S3 CCD, with a 200 kilosecond exposure time. In order to separate out 
potential calibration issues from the focus of this particular study, these response matrices 
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implement a constant response over the detector, and in generating the simulated X-ray 
images, the evolved zero-redshift clusters are shifted to a redshift sufficient to fit -R500 within 
the 16 arcminute field of view. Because of this step, though the clusters in the sample are 
distributed over a range of physical sizes, they all have approximately the same apparent 
size projected onto the image plane and into the simulated observations. 



2.3. M-T x Relation 



In this paper, cluster mass refers to M^nn, while Ty refers to our estimated "spectral- 
fit temperatures" obtained using XSPEC (jArnaudlll996l ) to fit single-temperature MEKAL 
plasma models to various energy bands in simulated X-MAS X-ray observations taken from 
the numerically-simulated clusters. The steps involved in the process are described in detail 
in §0 

Figure [3] shows the mass-temperature relation based on our sample of simulated clusters 
and estimated spectral fit temperatures, mea sured in a broad 0.7-7 keV band within an 
annular aperture from 0.15i?soo to -R2500 as in 



power- law form 



M = M 



Cavagnolo et all (120081 ). The best fits to the 
TX 1 (1) 



2keV, 

have the coefficients M = 1.09 ± 0.01 x 10 14 h~ l M & , a = 1.57 ± 0.06 for the full sample 
of clusters. For the subset of clusters with T 2 -7> 2.0kev, the best-fit parameters are 
M = 1.06 ± 0.02 x 10 14 hr x M Q , a = 1.71 ± 0.07. As is generally the case for simulated 
clusters, the power-law indices of the mass-tempe rature relations found here are co nsistent 
with cluster self-similarity and the virial theorem ( IKaiserlll986l ; iNavarro et al.lll995f ). These 
relationships have scatter, which we characterize by the standard deviation in log space a\ n u 
about the best-fit mass at fixed temperature T x . We find a\ n M = 0.10. 



3. Analysis 
3.1. Filtering 



We use X-MAS to simulate X-ray observations — provided as standard Chandra event 
files — then subject them to a series of r eduction steps using th e Chandra Interactive Analysis 
of Observations package (CIAO) v4.1 (IFruscione et al.ll2006l ). All of the event files have at 
least 100K counts, while those for the most luminous clusters have over half a million counts. 
Simulated observations of this quality provide us with an opportunity to address the question 
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Fig. 3. — Mass-temperature (M-Tx) relation for the clusters in our sample. Clusters with 
ksT < 2 keV are plotted in purple open triangles, while clusters with k B T > 2 keV are 
plotted in blue open squares, and the four massive D09 clusters in green open circles. Mean 
relations are plotted with solid lines for the whole sample, and dashed lines for clusters with 
ksT > 2keV. Finally, note that average temperatures are taken within a core-excised annulus 
whose outer radius is placed at i?25oo- [See the electronic edition of the Journal for a color 
version of this figure] 
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of how much of the scatter in our temperature-structure statistics is intrinsic. Consequently, 
our first reduction step filters each raw event file into a new set of event files by randomly 
sampling the events it records. The number of counts in each file is its "count level," and 
in order to cover the space of typical Chandra archival observations, we apply the filtering 
step to each of the raw event files four times, at count levels 15K, 30K, 60K, and 120K. 
These count levels roughly map to the range between observations of short duration or of 
relatively low surface-brightness objects to observations of long duration or of relatively high 
surface-brightness objects. 



3.2. Cool Lump Excision 

An example of a filtered event file appears in Figure H] as a surface-brightness image. 
It displays a common feature of numerical simulations of galaxy clusters, which is the pres- 
ence of relatively dense, cool, and metal-rich substructures that continuously undergo mass 
accreti on and have not yet come into thermal equilibrium with the hot ICM surrounding 



them ( IBorgani Kravtsovll2009f ). These bright point-like spots, or "cool lumps", typically 
are associated with the dense cool cores of smaller halos that have merged with the primary 
halo. Generally regarded as an unphysical artifact of numerical simulations, at least inso- 
far as their temperatures, densities, and co ncentrations are concerned, commonly they are 



excised before fur ther analysis is conducted (IRasia et al.ll2006l ; iPiffaretti fc Valdarninil 12008 



Naeai et all 120071 ) . 



In order to study the effect that excising cool lumps has on measures of temperature 
substructure, we produce from the originals several new event files whose cool lumps have 
been excised to an increasing degree. The CIAO wavdetect tool identifies peaks in the 
photon distribution by correlating an event file's image with a sequence of "Mexican-Hat" 
wavelet functions of differing scale sizes, measured in pixels, then generates a source list 
with associated region files. We use its default sequence of wavelet scale sizes in this analysis 
(2 and 4 pixels), though we apply the tool four times to each filtered event file, each time 
incrementing the multiplicative factor by which the source regions are scaled. In CIAO 
wavdetect, the parameter governing this multiplicative factor is ellsigma and in our study 
ranges between to 3. A value of for ellsigma is equivalent to "no masking" while a value 
of 3 corresponds to what we call "full masking". The effect is to produce versions of each 
cluster observation with a range of masking. Note that the wavelet scale size is a constant 
fraction of the cluster size because all the clusters have been redshifted so that .R500 fits in 
the field of view. 



We finish the extraction phase of our analysis with the following steps. First, for every 
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Fig. 4. — X-MAS simulated X-ray surface-brightness images for the same cluster whose 
spectrum is presented in Figures [1] and [2j The aperture is set at -R2500 and a core region of 
size 0.15i?25oo is excised. Note the surface-brightness peaks in the image. In this study these 
are detected automatically by the CIAO tool wavdetect, which generates source regions 
identified in this figure by red ellipses. Labeled "cool lumps" in this analysis, they are 
subjected to varying degrees of excision, from no masking (left) to full masking (right), as 
described in the text. [See the electronic edition of the Journal for a color version of this 
figure.} 
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excised file we apply an aperture of -R2500 and excise the central 0.15 -R2500 region, generating 
new copies o f the event files. This step, including the centering algorithm, matches the 
procedure in ICavagnolo et al.l (120081 ) so that regions measured in the simulated clusters 
correspond to those in the Chandra archival observations. Next, each fully-processed event 
file file is extracted into a standard pulse-invariant (PI) spectral file binned so as to have at 
least 25 counts in each energy channel. The extraction phase is complete when each original 
raw event file generates PI spectral files suitable for spectral fitting in XSPEC. 



3.3. Spectral Fitting 

The PI spectral files are then fed into XSPEC vl2.5.0 for spectral fitting. In order to 
form the temperature ratio Thbr two fits are performed for each spectrum. The first fit is over 
the hard band from 2.0(1 + z)~ x -l keV, for which the (1 + z)~ x factor exists in order to shift 
the 2 keV cutoff from the observer's frame to the cluster's rest-frame. The simulated clusters 
occupy a range of redshifts because larger clusters are translated to higher redshifts in order 
to fit i?5oo within 16 arc-minutes, when creating artificial observations. The second fit is over 
the broad band, including all energy channels in the spectrum from 0.7-7 keV. Every fit is 
made by minimizing the x 2 statistic for a single-temperature MEKAL model multiplied by a 
warm-absorber (WABS, to account for Galactic absorption). The Galactic column is fixed to 
Nh = 5 x 10 20 cm -2 and the metallicity to 0.3 Solar, leaving the temperature of the emission 
component, its H density (although this makes no difference), and its normalization as the 
only free parameters. 



3.4. Quantifying Temperature Structure 

We adopt two spe ctral measurements of temperature structure. The first is the tem- 



perature ratio Thbr of ICavagnolo et al.l (120081 ) and is found in the following way. For each 



simulated cluster in our sample, for each count level (15K, 30K, 60K, 120K), and for each 
cool lump masking level (0, 1, 2, 3), we find a spectral-fit temperature in the hard band and 
the broad band and form the temperature ratio Thbr- Emission from cooler, line-emitting 
metal-rich parts of the ICM ought to be excluded from the hard band, so we expect that a 
Thbr value greater than unity signals the presence of merging sub-clumps. 

Our second measure of temperature, which we call the "cool residual," (RES coo i) com- 
pares an observation's actual broadband count rate to a model-predicted count rate for a 
spectral model fit only over the hard band. Again, we determined a spectral-fit temperature 
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for each combination of cluster, count level, and masking level, except that in this case we 
performed only a hardband fit. Since emission from cooler components should be excluded 
from this band, in general we achieve good fits even when the count level is sufficiently high 
that broadband fits may formally have large reduced x 2 values. From this hardband model fit 
we estimate the corresponding broadband count rate and calculate its percentage deviation 
from the actual broadband count rate. While single-temperature systems will have actual 
count rates that are essentially the same as their model count rates, the introduction of 
cooler, luminous substructure components should generate an excess broadband count rate 
relative to the model. This measure of temperature substructure avoids the uncertainties 
associated with a single-temperature fit to the broad band. 



4. Results and Discussion 

Having calculated spectral-fit temperatures, T H br, and the cool residual for all of the 
clusters in our sample, we have the means to study how well temperature structure correlates 
with departures from the mean M-Tx relation, and how well these measures can be used to 
obtain better mass estimates. We will also examine how our results depend on the removal 
of the cool lumps. 



4.1. Thbr from Simulations. 



We present in Figure the temperature ratios T H br for our simulation sample, plotted 
as a function of t he br oadband temperature fit Tq.7_7. This figure is similar to Figure 8 in 



Cavagnolo et al.l (120081 ). which shows that the mean value of Thbr observed among clusters 
in the Chandra archive is (T H br) = 1-16, with a standard deviation of ct TH br = 0.10. We 
remind the reader of import ant differences be t ween the two samples being considered. The 
Chandra archive sample of ICavagnolo et al.l ( 120081 ) contains clusters most of which have 
ksT > 3 keV, whereas most of the clusters in our sample have ksT < 3 keV. Nevertheless, 
even with that caveat the temperature ratios of the simulat ed clusters are d i stribu ted in 
approximately the same way as are those in the real sample of ICavagnolo et al.l ( 120081 ). with 
(T H br) = 1.12 and (t TH br = 0.11, provided the cool lumps are not excised (circles in Figure 
|5]). When the cool lumps are excised (squares in Figure |5]), the mean and variance of Thbr 
diminish, with (Thbr) = 1-07 and othbr = 0.07. These values are for the full 120K count 
data, though the full set of (Thbr) for the four ellsigma values and four count values are 
tabulated in Table [TJ This table provides mean values for Thbr and its variance for all of 
the clusters in the sample, and for a subset whose T 2 .o_7 value is greater than 2 keV. Much 
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of the variance derives from the lower-temperature clusters, and as they are removed the 
variance in T H br drops significantly, especially when full masking is applied. 

Figure E] presents a similar effect for RES coo i. Here, we see that the excision of cool 
lumps again reduces the mean and variance of the temperature substructure measure, al- 
though in this case the effect is more dramatic. Evidently, the presumably unphysical cool 
lumps in simulated clusters may be necessary to reproduce the distribution of quantitative 
temperature substructure measures found in real clusters. This is a subject which we return 
to in Section 14.31 



Our original motivation for conducting this study was to determine if temperature struc- 
ture, as quantified by the temperature ratio T H br, correlates with and can be used to correct 
for departures from the mean mass-temperature relation M-Tx- In order to test this idea, 
we define the "mass offset" at fixed temperature by the relation 



where M is the cluster's actual mass, and M pre d is the mass predicted from the mean M-Tx 
relation. Figure [7] plots the mass offsets for our simulated clusters as a function of Thbr, for 
the case in which predicted masses are derived from the M-Tx relation for the broadband 
spectral fit temperature. Values of T H br calculated both with cool lumps excised (squares) 
and without excision (circles) appear in this figure. Error bars on Thbr are omitted for 
clarity. While there is some correlation, such that clusters with more temperature structure 
(larger T H br) tend to be more massive than predicted by the mean M-Tx relation, the trend 
is weak and has substantial scatter. Excising the cool lumps weakens the trend further. 
Figure [8] shows the same comparison for the RES coo i measure instead of Thbr, and in this 
case we again find that excising the cool lumps has an even larger effect for Thbr although in 
both cases accounting for temperature substructure does not greatly reduce scatter in mass 
offset. 



We now examine the decline in Thbr and its variance as cool lumps are excised, begin- 
ning with Figure [91 which focuses on Thbr- Here, we plot the statistics standard deviation 
0"thbr as a function of the two dimensions along which we adjust our analysis pipeline, with 



4.2. Temperature Structure and Scaling Relations 
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4.3. Effects of Masking Strategy 
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Fig. 5. — Thbr plotted against T 7 _ 7 for simulated X-MAS observations that have approxi- 
mately 120K counts. Hydrodynamic simulations may produce spurious over-condensations of 
cool gas. These cool lumps can be excised using the CI AO tool wavdetect, whose aggressive- 
ness can be controlled via its ellsigma parameter. Squares correspond to single-temperature 
MEKAL fits whose underlying observations are processed by wavdetect with ellsigma set to 
3, and circles correspond to those with ellsigma set to 0. Finally, note that the 4 massive D09 
clusters are denoted by filled symbols rather than by open symbols. The solid line represents 
the mean when cool lumps are not removed, while the dashed line represents the mean when 
they are removed. The dotted line indicates Thbr = 1- [See the electronic edition of the 
Journal for a color version of this figure.] 



- 16 - 



120 





100 - 




80 - 


o 




c 


60 - 


w 




53 


40 - 




20 - 




10 



k B T 2.0-7.0 ( keV ) 



Fig. 6. — RES cooi plotted against T .7_ 7 for simulated X-MAS observations that have ap- 
proximately 120K counts. This statistic tracks the percent excess count rate of the actual 
observation with respect to count rate for a MEKAL model fit to just the hard band. Squares 
correspond to single-temperature MEKAL fits whose underlying observations are processed 
by wavdetect with ellsigma set to 3, and circles correspond to those with ellsigma set to 0. 
As in the previous plot, the 4 massive D09 clusters are denoted by filled symbols rather than 
by open symbols. The solid line represents the mean when cool lumps are not removed, the 
dashed line represents the mean when they are removed, and the dotted line indicates zero 
residual. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 7. — Relationship between Thbr and mass offset S In M(Tx) from the mean M-Tx 
relationship for the clusters in combined B04+D09 simulation sample. The temperatures in 
this relation are our broadband spectral-fit temperatures. Squares correspond to simulated 
X-MAS observations that are processed by wavdetect with ellsigma set to 3 (full masking), 
while circles correspond to observations that have ellsigma set to 0. [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 8. — Relationship between RES coo i and mass offset S In M(Tx) from the mean M-Tx 
relationship for the clusters in combined B04+D09 simulation sample. The temperatures in 
this relation are our broadband spectral-fit temperatures. Squares correspond to simulated 
X-MAS observations that are processed by wavdetect with ellsigma set to 3, while circles 
correspond to observations that have ellsigma set to 0. [See the electronic edition of the 
Journal for a color version of this figure.} 
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the top panel devoted to the masking strategy, and the bottom panel devoted to the count 
level. Focusing on the top panel, we see that increasing the ellsigma parameter of the 
CIAO wavdetect tool from to 3 reduces 0thbr from 0.15 down to 0.13, when the observa- 
tions are relatively "poor" (with ~15K counts). With higher-quality observations of 60K or 
120K counts, the decline is in othbr greater, going from 0.12 down to 0.07 as the ellsigma 
parameter rises from to 3, and ctthbr for full ma s king e nds up being significantly less than 



the value of 0.10 observed in the ICavagnolo et al.l (120081 ) sample. 



Figure [TO] helps show why the temperature ratio T H br and its variance decline as the 
cool lumps are removed. It plots the sample average of the relative change in temperature 
as the CIAO wavdetect ellsigma parameter is increased. As more of each cool lump is 
excised, both the broadband spectral fit temperature, and the hardband fit temperature 
increase. However, the increase is significantly larger for the broadband temperature, as the 
cool lumps' contribution to the flux is already largely excluded from the hardband fits. The 
top panel in this figure shows this effect for spectra with approximately 120K counts, while 
the bottom panel is for spectra with approximately 15K counts. 

Another way of visualizing the effect of more aggressive masking is depicted in Figure 
rTTi The top panel in this figure is similar to Figure |6] in that it occupies the To.7_7-Thbr 
plane, with Thbr plotted as a function of T .7_7 for our simulated clusters. Arrows illustrate 
the shift in Thbr and in T .7_7, as the cool lumps are excised. Some of the clusters experience 
very large shifts in both quantities, while others experience no shifts at all. The former are 
associated with clusters that have many well-defined and bright cool lumps, while the latter 
correspond to those clusters that are completely free of cool lumps. The bottom panel in 
this figure shows just the size of this shift for the Thbr statistic, from which we can see that 
some clusters indeed have a shift of precisely 0. Again, these are clusters whose simulated 
X-ray observations are unchanged after applying the CIAO wavdetect tool because it finds 
no sources to mask out. 



4.4. Substructure Measures and T H br 



Various researchers have established a cl ear correlation between morphological mea- 



sures of cluster dynamical state and 8 In M (I Jeltema et al.l 120081 ; IVentimiglia et al.l 12008 



Piffaretti fc Valdarninil l2008t lYang et al.l 120091 ) , a corr elation which is less ap parent in our 
study of spectral measures of substructure. Similiarly, ICavagnolo et al.l (120081 ) found a cor- 
relation between the temperature ratio and structure for their Chandra sample, in that 
merging events are associated with elevated Thbr- To probe this issue further we focus on 
Thbr and compare it to several metrics for cluster substructure. These are the centroid 
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Fig. 9. — {top) Decline in the Thbr standard deviation <7thbr as wavdetect's ellsigma pa- 
rameter ranges from to 3. These are plotted for a family of simulated observations of 
increasingly higher quality, from approximately 15k counts to approximately 120k counts. 
{bottom) Decline in othbr as simulated X-MAS observations go from lowest-quality (approx- 
imately 15k counts) to highest- quality (approximately 120k counts). These are plotted for 
a family of simulated observations with wavdetect's ellsigma parameter ranges from to 3. 
[See the electronic edition of the Journal for a color version of this figure] 
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Fig. 10. — Relative change in spectral fit temperatures as the cool lumps are removed 
from the simulated clusters for high-quality simulated observations, (top) Observations with 
approximately 120K counts, (bottom) Observations with approximately 15K counts. Note 
that these are for observations in an aperture corresponding to R 2 5oo-[See the electronic 
edition of the Journal for a color version of this figure.} 
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Fig. 11. — ( Top) Shift in the T .7_7-Thbr plane as wavdetect's ellsigma parameter ranges from 
to 3. Clusters from the B04 sample are depicted with solid arrows, while clusters from the 
D09 sample are depicted with long-dashed arrows, broadband spectral-fit temperatures are 
estimated within an aperture set at i?25oo- Notice that the distribution of most of the points 
shifts significantly to higher spectral fit temperatures and smaller temperature ratios as cool 
lumps are excised. Notice also that some simulated clusters are unafflicted by cool lumps, 
so that their net shift is 0. (Bottom) The Thbr component (y-axis) of the shift presented in 
the top panel. Note that these are for simulated X-MAS observations of maximum quality, 
having approximately 120K counts. [See the electronic edition of the Journal for a color 
version of this figure.} 



-23 - 



variation w, the axial ratio 77, and the power ratios P20 and P30 (see IVentimiglia et al.ll2008 
and references therein). The centroid variation w measures the skewness of a cluster's two- 
dimensional photon distribution by calculating the variance in a series of isophotes for the 
cluster surface-brightness map. The axial ratio rj measures a cluster's elongation, which 
tends to increase during merger events. The power ratios P20 and P30 decompose a cluster's 
surface-brightness image into two-dimensional multipole expansions, capturing different as- 
pects of a cluster's geometry. P 20 relates to the ellipticity in an image and is similar to the 
axial ratio 77, while P30 measures the "triangularity" in an image. 



In IVentimiglia et al.l (120081 ) we calculated these morphological substructure metrics for 
a superset of the B04 sample of simulated galaxy clusters and compared them to departures 
from the M-T x relation (see Fig. 8 in that paper). Here we use the same substructure 
measures for the 114 B04 clusters for which we are able to calculate Thbr, and co mpare the 



results. These are presented in Figure [12j Whereas in IVentimiglia et al.l ( 120081 ) there is a 
clear correlation between morphological substructure and S In M, we find in this study that 
there is little or no correlation between substructure and Thbr for our B04 sample. 

In order to understand how temperature structure and morphological structure can be 
uncorrelated, we looked at four simulated clusters having approximately the same temper- 
ature (ksT ~ 3.2A;ey) and occupying the relative extremes in Thbr and in the centroid 
variation w. Two were selected for relatively low Thbr (— 1-1) but large variation in w 
(~ 0.02-1.0). Two were selected for large Thbr (— 1-4) but again large variation in the 
morphological structure parameter u (~ 0.1-1.0). These clusters are presented in Figure ITB1 
with surface-brightness contours overp lotted and with the associated values for Thbr and w. 
The two clusters in the left column appear more symmetric in their contours, while the two 
in the right column exhibit noticeable centroid shift. However, the morphologically apparent 
structure in the lower right cluster is not observed spectrally in the temperature ratio T H br- 
Evidently, neither morphology nor Thbr is a perfect measure of relaxation in our simulation 
sample, as it contains clusters with small centroid shift u and obvious substructure, as in the 
upper left of this figure. And, it contains clusters with large w that are nearly isothermal, 
as in the lower right of this figure. 



5. Summary 

We used a sample of galaxy clusters simulated with radiative cooling and supernova 
feedback, along with simulated Chandra X-ray observations of these clusters, to study tem- 
perature inhomogeneity as a signature of cluster dynamical state. Specifically, we adopted 
two methods of quantifying temperature inhomogeneity spectroscopically, the temperature 
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ratio Thbr ( IMathiesen h Evrardll200ll ; ICavagnolo et al.ll2008l ) and the cool residual RES coo i. 
The former is the ratio of a hardband X-ray spectral-fit temperature to a broadband tem- 
perature, and becomes greater than 1 for clusters whose ICM contains cool, over-luminous 
sub-components. The latter is the excess broadband count rate relative to the count rate pre- 
dicted by a model fit to the hard X-ray band. Though our simulated clusters are typically less 
massiv e and have lower temperatures than the Chandra archive clusters in ICavagnolo et al. 
(120081 ). we find that their temperature ratios T H br occupy generally the same distribution 
as the observed clusters. 

We also looked for an opportunity to combine Thbr and the cool residual with the mean 
mass-temperature relation to obtain better mass estimates than are achieved just with the 
scaling relation alone. We find, however, that while both T H br and the RES coo i are correlated 
with offset from the M-Tx relation, these correlations are weak, at least for this sample. We 
conclude that these measures of temperature inhomogeneity are not very effective at reducing 
scatter in the mass-temperature relation. 

Finally, we note a particular difficulty that arises when trying to use clusters from hy- 
drodynamic simulations to calibrate scatter-correction observables based on temperature in- 
homogeneity, such as Thbr- Historically, simulated clusters have tended to exhibit their own 
kind of "over-cooling problem" , in which dense lumps of luminous, cool gas associated with 
merged sub-halos appear. These cool lumps are often excised from simulated clusters before 
their global properties are measured, but masking them appears to make the remaining tem- 
perature structure of the simulated clusters overly homogeneous. This finding suggests that 
real clusters may have cooler sub- components, that are more diffuse and less concentrated 
than in their simulated count erparts. Some physical process, perhaps thermal conduction, 
turbulent heat transport (e.g.. iDennis fc Chandranll2005r iParrish. Quataert fc Sharmall2010l ; 
Ruszkowski fc Ohll201l[ ). or a more aggressive form of feedback, prevents cool lumps from 
forming in real clusters and might not completely eliminate those temperature inhomo- 
geneities. Newer simulations incorporate treatments of conduction and AGN feedback, as 
well as more accurate treatments of mixing, and it will be interesting to revisit these temper- 
ature inhomogeneity measures in simulated clusters when large samples of such simulated 
clusters become available. 
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Fig. 12. — Relationship between Thbr and four measures of substructure: the axial ratio rj, 
the centroid variation w, and the power ratios P20 and P20. Squares correspond to simulated 
X-MAS observations that are processed by wavdetect with ellsigma set to 3 (full masking), 
while circles correspond to observations that have ellsigma set to 0. [See the electronic edition 
of the Journal for a color version of this figure] 
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Fig. 13. — Surface-brightness contours for four clusters drawn from our simulation sample, 
illustrating that temperature ratio and centroid shift are both imperfect measures of relax- 
ation. All four clusters are at nearly the same temperature, with k B T ~ 3.2. The two in 
the left column exhibit little centroid shift w, while the two in the right column have cen- 
troid shift near the maximum for the sample. The two clusters in the bottom row have low 
Tubr while the two in the top row have higher T H br suggesting the presence of multiple 
temperature components. 



